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The generation of photons in a three dimensional rectangular cavity with two moving boundaries 
is studied by using the Multiple Scale Analysis (MSA). It is shown that number of photons are 
enhanced for the cavity whose walls oscillate symmetrically with respect to the center of the cavity. 
The non- stationary Casimir effect is also discussed for the cavity which oscillates as a whole. 
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I. INTRODUCTION 
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CN , When two perfectly conducting plates are placed close to each other, the attractive force appears between the 
parallel conductors due to the vacuum fluctuations as predicted long time ago by Casimir The corresponding 
\ vacuum energies and forces were static. Now let us assume that the right boundary depends on time. In this case, the 
[t I length of the cavity changes in time, L = L(t). The most evident manifestation of dynamic behavior is the dependence 
of the force on time. The modified force depends on time [2]. 



(N 



2AL{t) ■ 

More strikingly, when the right plate is in motion, it was theoretically predicted that photons are generated in the 
i' empty cavity, because of the instability of the vacuum state of the electromagnetic field in the presence of time- 
D ', dependent boundary conditions [1, 0, d, 0, 0, H, d, [13, [HI, ■ A number of virtual photons from the vacuum are 
converted into real photons. This phenomenon is known as dynamical Casimir effect or motion-induced radiation. 
\ However, there has been no experimental verification for this effect up to now because of the technical difficulties. 
^— I . There are a few proposed experimental for the detection of photons [ij, Uj, . As it was discussed in the literature, 
K>" ■ the best way to observe this effect is to vibrate one of the wall with one of the resonant field frequencies. A one 
dimensional cavity with two perfectly parallel reflecting walls, one of which is motionless and the other oscillating 
with a mechanical frequency equal to a multiple of the fundamental optical resonance frequency of the static cavity, 
, has been used as a simple model to study the dynamical Casimir effect [1, S 0, • 

r4 ■ Li(t) = 0, L2(t) =L(l + e sinf7t) , (2) 

o ; 

00 , where the constant is the external frequency, e is a small parameter, the constant L is the initial length of the cavity 

■ and Li{t) and L2{t) are the positions of the right and the left walls at time t, respectively. The cavity is motionless 

^ ' initially and that at some instant one mirror starts to oscillate resonantly with a tiny amplitude. 

• ^ . Calculating the number of generated photons is a difficult task since one has to solve the wave equation with the 

/\ time-dependent boundary conditions. A lot of technique have been developed to deal with the problem. For example, 

^ , averaging over fast oscillations [l^, , multiple scale analysis [l^ , the rotating wave approximation [l9| , numerical 

■ - - ' techniques [20| are applied to the dynamical Casimir problem. 



The case of cavities with two moving mirrors was studied by few authors [21|, [2^, l23|, |2J, |25|, |2a, l27|, |28| . Compared 
to the situation with a single oscillating mirror, it was found that radiation is resonantly enhanced when the cavity 
oscillates as a whole, with its mechanical length kept constant and when the cavity oscillates symmetrically with 
respect to the cavity center. The radiation emitted by two oscillating walls in one dimension was studied by Dalvit 
and Mazzitelli using renormalization group method [23], by Dodonov using the method known in the theory of 
parametrically excited systems [25|] and by Lambrecht et al. using the scattering approach [28| . 

In this paper, we will study two moving boundary problems for a rectangular cavity resonator in 3 — D. We will use 
multiple scale analysis (MSA) which provides us with a solution valid for a period of time longer than that of the 
perturbative case, which is not suitable for this problem since it breaks down after small time because of the resonance 
terms. We will investigate the two cases. In the first configuration, we will consider a rectangular cavity resonator 
whose right and left walls in the x direction move in exactly the same way. 

Li{t) = e LsinQt , L2{t) = L (1 + e sinQt) , (3) 

Initially, the length of the oscillating cavity is given by L2{t) — Li{t) = L. As time goes on, the length of the cavity 
is not changed. This type of motion corresponds to the cavity oscillating as a whole. 
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Secondly, we will consider a rectangular cavity resonator whose right and left walls in the x direction move opposite 
to each other. 

Li{t) ^ -e LsinQt , L2{t) ^ L (1 + e sinQt) , (4) 

We will show that the number of generated photons are enhanced for this system with the parametric resonance case, 
in which the frequency of the wall is twice the frequency of some unperturbed mode, say = 2ujk. 
This paper is organized as follows. Section II studies the field quantization in the case of moving boundaries in one 
dimension. Section III reviews and applies the multiple scale analysis. Section IV studies the dynamical Casimir 
effect for a three dimensional cavity. Section V considers the cavity whose walls are oscillating symmetrically with 
respect to the center of the cavity. Finally the last section discusses the enhancement of generated photon number . 



II. FIELD QUANTIZATION WITH FIXED LENGTH 

Consider a one dimensional cavity formed by two perfect conductors. The right and the left walls oscillate according 
to the formula given by ([2|). The cavity oscillates as a whole. The field operator in the Heisenberg representation 
^{x,t) obeys the wave equation (c = 1) 



The boundary conditions are given by 



$(Li,i) = $(L2,i) = , (6) 



which describes the moving boundary problem in field theory. The problem looks simple, since the wave equation 
is the same as whether the boundaries are moving or not. However, the moving boundary conditions render the 
equation unsolvable by the usual means. Although the solution of the wave equation is easy to be found and well 
known, finding the exact solution of the same problem endowed with the time-dependent boundary conditions is very 
difficult and not known except for some special cases. 

To solve the problem, we will transform the moving boundary conditions to the fixed boundary conditions. Let us 
introduce a coordinate transformation as 

m ^ . (7) 

As a result, the new time-independent boundary conditions for ^{q,t) read 

<i>(g-0,i) = $(g=l,t)=0. (8) 

With this coordinate transformation, the wave equation ([5]) is changed. Let us find the new form of the wave equation. 
Under the coordinate transformation ([7]), the derivative operators transform as 



91 - ^ dl ■ (9) 
d d 

Here we use the notations 9t = — , 9„ = — . By using the explicit forms of Li(t) and L2(t) the transformation 

dt dq 

of the time derivative operator up to the first order of e can be approximated as 

^df ~e{2fl cosflt dtdq - V? sinl7i 9, ) . (10) 
Substituting these into the wave equation, we get 
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In the right hand side of the equation, the term in the parenthesis is due to the time dependent boundary conditions. 
Now, we don't have to deal with the time dependent boundary conditions. However, we are left with a new equation. 



3 



Let us now solve the above equation subject to the fixed boundary conditions ([8]). The field operator can be expanded 
as 

Hq,t)^J2 {bk^kiq,t)+bin{q,t)) , (12) 
k 

where 6[.and bk are the creation and the annihilation operators, respectively and 'i'k{q,t) is the corresponding mode 
function. We will follow the approach given in [2^, [30, [3l| to find the explicit form of functions \E'fc(x,t) . For an 
arbitrary moment of time, the mode function satisfying the boundary conditions is expanded as 

*fe(g,t > 0) =J2a';Xt) sin(n7rg) , (13) 

n 

Let us substitute it into the equation (|lip in order to find a^(t). Then, multiply the resulting equation with sin(rn7rg) 
and integrate over q from to 1. If wc use the orthogonality relations, we get an infinite set of coupled differential 
equations for a'^ after some algebra. 

+ uj^^ a™ = e cos ilt ^ g„„ dj; - 2^^ sin Qt ^ g,,^ (14) 

mn(l - (-1)™+") 

where ujm — mi: IL and the antisymmetric coefficient is given by gnm = 5 for m ^ n. 

In the next section, we will solve the above equation. To do this, we prefer to use multiple scale analysis method. 



III. MULTIPLE SCALE ANALYSIS (MSA) 



Conventional weak-coupling perturbation theory suffers from problems that arise from resonant terms in the per- 
turbation series. The effects of the resonance could be insignificant on short time scales but become important on 
long time scales. Perturbation methods generally break down after small time whenever there is a resonance that lead 
to what are called secular terms. In the equation (|14p . this happens for those particular values of external frequency 
f2 such that there is a resonant coupling with the eigenfrequencies of the static cavity. To avoid such problems, 
we will use Multiple-scale analysis (MSA), a powerful and sophisticated perturbative method valid for longer times. 
Multiple-scale perturbation theory provides a good description of our system. 

The trick is to introduce a new variable r = et. This variable is called the slow time because it does not become 
significant in the small time. The functional dependence of aj^ on t and e is not disjoint because it depends on the 
combination of et as well as on the individual t and e. The time variables t and r are treated independently in MSA. 
Thus, in place of a^(t), we write a%i{t^ et). Let us expand af„ in the form of a power series in e 

alit) = af„(°'(t,r) +6 aJ^W(t,r) -He^ a,T'(i,r) + ... . (15) 

To this end, we change the independent variable in the original equation from t to t. Using the chain rule, we have 
± d_ d_ 

Up to the first order of e, the derivatives with respect to the time scale t are given by 

hi = dl a^(") + e (2drdt a^^^ + d^, a^^^) , (16) 

where dot denotes time derivation with respect to t as usual. Let us substitute these into the equation (|14p . Then, we 
see that our original ordinary differential equation is replaced by a partial differential equation. It may appear that 
the problem has been complicated. But, as will be seen below, there are many advantages of this method. To zeroth 
order in e, we get a well-known equation in physics. 

dt(°)-Hc.^aJ^(°)=0. (17) 

To first order in e, we obtain the following equation. 

afi^) + c^^ajj(i) = -2dtdr af;(°) + ^ cos{nt) ^ g„„ d^(°) - sm{nt) ^ g„™ a^(°) (18) 
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The solution of the former one can be found easily 

a';i°Ht,T)^AUr) e-^--'+Btir) e^'"* . (19) 

Note that A^{t) and B^{t) are not constants but functions of the slow time scales r. The initial conditions are given 
by 

Btir = 0) = (20) 

Let us now solve the equation p^ . We look for the oscillations of the cavity that could enhance the number of 
generated photons by means of resonance effects for some specific external frequencies. To enhance the number of 
photons, let us now assume the resonance condition, = p n/ L, where p = 1, 2, .... It is well known that, in the 
resonance conditions, the number of generated photons grows very much in time. 

Let us substitute the zeroth order solution (flQ)) into the right hand side of the equation (jlSp and then use the following 
relations: 2i sin fit = (e'^* — e^*^^*), 2 cos fit = (e*^* + e~*^*). It can be seen that the right hand side contains terms 
that produce secular terms. For a uniform expansion, these secular terms must vanish. In other words, any term with 
g±iuj,„t ^YiQ right-hand side must vanish. If not, these terms would be in resonance with the left-hand side term 
and secularities would appear. After imposing the requirement that no term e^"^'"* appear, we get 

m ' ^p+m.va p+ra ^rn—p,m m~p ^p~rn.m p—m ^ V / 



where Gtj is defined as 

In the similar way, the fact that no secularities should arise from the term with e'""* leads to 



G^, = ^ n . (22) 



T m p+m,m p+m ^ rn— p,m rn—p p— m,m p—rii v 

To this end, let us give the formula for the number of generated photons 0, H^l 

(iV„>^5]2c.„|B^r (24) 

k 

We will now analyze the solutions of the equations (|21l23p for a given p. 



A. Analysis of Solution 

As a special case, let us study the above equations when m — p. In this case, the equations (|21l23p are reduced to 
the following simple ones. 

9rA'p + G2p,p Ap = (25) 

dr Bl + ,p B% = (26) 

To find the solution, we should also find the corresponding differential equations for A2p,J52p- Substituting m = 2p 
in (|21I23P gives 

9t B^p + G^p p B^p — Gp 2p Bp — (28) 

Note that since p — m^p — 2p<0, the term with Bp_^ vanishes. As can be seen easily, we need the equations 
A\ and i?| to solve these differential equations. In fact, there is no upper cutoff. However, these equations are 
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not coupled equations. A'^ and are not coupled to each other. So, we say that no photons with ujp,uj2p, ■■■ are 
generated for a given external frequency Q = irp/L. But this doesn't mean that photons are not generated. For 
example, photons with ujp+i are produced. To calculate the number of photons in this mode, we should solve (|21l23p 
when m = p + 1. As was seen above, we are left with infinitely many equations since there is no cutoff. This is due to 
the fact that the spectrum of a one-dimensional cavity is equidistant. Intermode coupling produces resonant creation 
in the other modes. So, the spectrum docs not have an upper frequency cutoff. So, multiple scale analysis does not 
estimate the number of generated photons in one dimension. 

To sum up, photons are produced resonantly in all modes except m — p, 2p, .... In the next section, we will study the 
dynamical Casimir effect for a three dimensional cavity whose spectrum is not equidistant. 



IV. DYNAMICAL CASIMIR EFFECT IN 3-D 



So far, we have restricted ourselves to the one dimensional case. We will now study the dynamical Casimir effect for 
the three dimensional geometries. We will show that multiple scale analysis works very well to calculate the number 
of generated photons. 

Let us firstly define our problem. Consider a rectangular cavity resonator with perfectly conducting walls. Initially 
the three of the them are placed at x = 0, y = and z — while the other three walls are at x — Lx,y — Ly and 
z = Lz- At time t = Q, the cavity starts oscillating in the x-direction as a whole. The positions of the six walls at any 
time are given by 

Li^t) = e L^smnt , La^W = (f + e sin f^O , Liy(t) = , L2y{t) = Ly , Li^t) = , L2z{t) = . (29) 

Here the constant f2 is the frequency of the oscillation and e is a small parameter. The volume of the cavity resonator 
is constant in time and given by = LxLyLz. 

As in the one dimensional case, we will study with the coordinates transformed from the fixed ones to the moving 

ones, (x, y, z) (g, y, z). The relation between q and x is given by g = (x — LiJ)/Lx. 

We will firstly study the dynamical Casimir effect for the scalar field and then for the vector field. 



A. Scalar Field in 3-D 



In this section, we will apply MSA to the scalar field in three dimensions. The cavity resonator oscillates in the 
x-direction ((29)) . So, the scalar field operator subject to the following boundary conditions 



y, z, t) = ^{L2x,y, z, t) = $(x, 0, z, t) = $(x, Ly, z, t) = $(x, y, 0, t) ^ $(x, y, L^, t) = , (30) 

which describes the moving boundary problem in three dimensions. The field operator in the Heisenberg representation 
^^(x, y, z, t) obeys the wave equation (c = I) in three dimensions 

We can rewrite this equation in the moving coordinate systems, {q,y,z). The transformation of the time derivative 
operator is given in pop . Hence, the wave equation is transformed to 

In one dimension, the expansion of the solution of the wave equation was given by the formula (|I2p. In three dimension, 
it has obvious generalization. The scalar field operator can be expanded as 

^q,y,z,t)= J2 bk^kyk^ 'i'k^kykA<l,y,z,t) + H.C. . (33) 

where bk^k k^ is the annihilation operator and VP/c^fe k^iQTy, z,t) is the corresponding mode function. Let us assume 
that, for an arbitrary moment of time, the explicit form of function ^'^^fc^fc^ (q, y, z, <) is expanded as 

*fc.fe„fc.('Z,y,^,^>0)= J2 "rtn/iW sin(n,^<z) sin(^y) sin(^z) . (34) 
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Substituting these into the transformed wave equation ((32|) and using the orthogonaHty relations, we get an infinite 
set of coupled differential equations for a^^^j"^^^(i) after some algebra. 



■■k^kyk^ I 2 



yna'^^^jX = ^ I 4r!cos(m) 5n.™. ailif^l - 217^ sin(m) ^ g„^™^ a^^^^^^ ) (35) 



where w™^„ „^ = i^'^i ml/Ll + nl/Ll + nl/Ll ) and .g„^™^ = 2 2 ^ 



equation is the three dimensional generalization of (jl4p . 

We will apply the multiple scale analysis to solve this equation. Fortunately, no need to go into the detail. The 
generalization of the equation (jlSp to three dimensions is straightforward. 



fjk^kyk^ii) , 2 k^kyk,{i) _ -28 dt n'^^^y'^''-^^ -\- 4Qcn-^(Qf) o 

-2n'sm{nt) ^ gn.n.. atltlifK (36) 



where 



a'^^riyu^!'^ = A'^^riyk (r) exp (-z c^™,„,„, t) + B W;.^ (^) exp {i w„^„^„^ i) , (37) 
^m'^n'^n (''') ^-i^d ^m'^n'^n (''') depend ou the slow time parameter r. The initial conditions are given by 

31 y z X y z 

^tny'^nAr = 0) = ^ 4x,nx 4„,n, ; 

5':n"„t(^-0) = 0. (38) 

We will derive an equation like (|21l23p . Substituting ([ST]) into and eliminating the secular terms from the 
equation, we get 

'^T m^riyn^ (nj^nyn^),(mxnyn^) n'^nyU^ {n'^riyn^) ,(m^nyn^) n'^nyTi^ (ii'J.'nyn^),(m^nyn^) n'^'riyn^ " 

n rfkxkykz r^k^kyk^ rtk^kyk- Akj^kyk^ ^ (A{]\ 

^r-i^rn^nynz (n'^nyiiz) .(m^nyn^) n'^nyn^ ' (n'^nyn^) Xm^Uyn^) n'^riyU, ' (n'J,' nyii^) Xm^nyn^) n'^'nyn^ ~ " I 

where the three dimensional generalization of the definitions ([22|) is given by 

^ "F ^WrxTiyrij „ , , 

Gj . , . — il Qr m ■ 41 

(r^,ny,nz),(m^,ny,nz) 9, , i"x"'x V y 

It should be noted that n'^ , n" , n"' are positive integers and they satisfy 

^n'^TLyTiz — ^" ^7n^ny7iz 7 

^n'^TLyTiz ^ ^ ^nixTiyrLz ; 

^n'^'riyUz ^ ^m^nyUz ■ (4^) 

In one dimensional case, ^771, ^m^pi ^7n=p2p; ■•■ and B^^^ -^m^p? -^7n=p2p; ■■■ are strongly coupled to each other. However, 
in three dimensions, only a few modes are coupled to each other. This is because there are only a few positive integers 
n!^,n",n^" satisfied by the equations So, it is possible to solve the equations p9l40p exactly since only a few 

modes are coupled. 

In what follows, we will give some specific examples. 



1. Examples 



For the simplicity, assume that the cavity is cubic, = Ly = Lz- We have a freedom to choose the coupled 
modes. For example, assume that inter-mode coupling occurs between (n"',ny,nz) and {mx,ny,nz). So we can 
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determine Q. from flS]) . As an example, we are interested in the following two modes: (1, 1, 1) and (2, 1, 1). 
n = (\/3 + V&)n/Lx- Let us solve the equations p9l40p for these modes. Hence, 

t^T^211 ""(14,1), (2, 14) -"111 ~ " ' 

drykj^kykz . ik^kykz ^ 

T-°211 + '-'(144), (244) ^111 ~ ^ ■ 

To solve the above differential equations, we also need the equations for drAi^i^i and 9i--Si,i,i- Using again 
we get 



Ct-^III 



G 



B. 



(2, 1,1), (1,1,1) -"211 



^a; ky 



0; 



T"-°lll + '-'(2, 1,1), (1,1,1) ^211 ~ " • 

The solution of the equations (|43l44p with the boundary conditions ((38l) can readily be found. 



Choose 

(43) 
(|39]40)1 . 

(44) 



111 
11 




4211 

^211 / \ 
where the constant A is defined as = G 



r211 

-c'lii 



Rill 
-^^211 



/ sinh(Ar) \ 



V 



'(1,1,1), (2, 1,1) 



(45) 



AV2ti^iii 



iinh(AT) J 



(2,1,1), (1,1,1) "(1,1,1), (2, 1,1) 

we find the number of generated photons for each mode 

(iVi,i,i) = sinh2(AT) 
(iV2,i,i) = sinh2(AT). 



^ With the help of number operator ([24] 



(46) 



(47) 



B. Vector Field 



Here, we will study the dynamical Casimir effect for the vector field in three dimensions. Maxwell's equations de- 
scribe electromagnetic waves as having two components, the electric field, E(x, y, z), and the magnetic field, H{x, y, z). 
Modes in a cavity resonator are said to be transverse. It is convenient to classify the fields as transverse magnetic 
[TM) or transverse electric {TE) according to whether E or H was transverse to the direction of oscillation. 
We will study the dynamical Casimir effect for the transverse electric modes. Before applying our formalism, let us 
first write TE modes in static case. In the mode TE, for t < the cavity is static, and each mode is given by 



A, = 

Ay ^ E exp (-iLj„^„^„^i) sin 



' cos [ — y 



Lx J \ L 



Tlx ,ny ,nz 



Az = 2^ Aq^ exp {-iuj„^nynj) sm \^ — — j j cos j , (48) 



where nx,ny,nz are positive integers. The constants Aq^ and Aq^ satisfy the coulomb gauge condition, 
Ao^Uy/Ly + Ao^Uz/Lz = 0. 

At time t — 0, the rectangular cavity resonator starts to oscillate in the x-direction as a whole. The positions of 
the six walls are given by ([29|) . Let us find the components of the field operator at any time. We will study with 
the moving coordinate systems (q, y, z). The field operator A(q, y, z, t) associated with a vector potential satisfies the 
transformed three-dimensional wave equation (c— 1) 

1 d'^A d^A d^A d^A / ^ ^ d^A ^2 . ^ dA\ 
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When t > 0, the solution of the components of the field operator may be expanded in terms of the orthogonal basis 
functions. 



A,. = 



Ay ^ a^i'ttr, {t) sin {n^nq) cos 



n-j. ,ny 



A.= Y. it) sin {n.nq) sin hfL) cos (^) 



(50) 



where the time dependent function a'^^^^^^ (t) is to be determined later. Let us substitute the equation ([50l) into the 
transformed wave equation Using the orthogonality relations, we obtain the dynamical equations. 



2 feajfeyfez 



e I 4fl cos(rJt) ^ 



'kx ky k z 



ky 



Sin(m) ^ g„^™^ "nxnyn. 



(51) 



As can be seen, this equation and the equation ((35| for the scalar field are the same. So, for the dynamical Casimir 
effect problem, both scalar and vector field cases can be treated in the similar way. Multiple scale analysis has already 
been applied to the scalar field. So, the equations (|39l40p are also valid for TE modes. We will study the photon 
generation for the TE modes by giving an example. 

We demand that the generated photons have the modes (2,2,1) and (3,2,1). These two modes are coupled if we 
choose r2 = (3 + VTa)^/ L^- 

After do the same things with the help of MSA we obtain. 



A kykz ) 

2,2,1 
A (kx kykz ) 

3,2,1 



-^2ci;2,2,l 

•y/2li^3,2,l 



4^,2COsh(A'T) 

4^^3Cosh(A'T) 



rjikxkykz) 

Tjikx kykz ) 
3,2,1 



'^[:-^;^'''^''''^' 5fe,,2sinh(A'r) 
^g.^.^).(3,^.i) 4^ 3sinh(A'r) 



(52) 



where A' = G^2i 221 ^221 321- Then, we calculate the number of generated photons ([24l 



(iV2,2,l) 
(^^3,2,l) 



sinh^(A'r) 
sinh^(A'T) 



(53) 



Here, we have performed analytical calculations to find the number of generated photons by using MSA. To this end, 
it should be mentioned that Ruser found perfect agreement between the numerical results and analytical predictions 
obtained by MSA (s^- 



V. ENHANCEMENT OF PHOTON NUMBERS 



So far, we have considered the cavity resonator oscillated as a whole. We will now study the case of symmetric 
oscillation with respect to the center of the cavity in the x direction. 

iix(i) = -e ixsinf^i , 12^1) ^ L^{1 + e BinVlt ) , Liy{t) = , L2y{t) = Ly , Li^{t) = , L2z{t) ^ L, . {54) 

In this case, the volume of the cavity changes in time. The two walls in the x direction move opposite to each other. 
For this configuration, only the scalar field will be treated. This is because, as was pointed above, the equations for 
'^ri'^'n (0 ^I'S the same for the scalar field and the vector field. In other words, the number of produced TE-mode 
photons equals the number of produced scalar particles in a three dimensional cavity. So, it is enough to study the 
dynamical Casimir effect for the scalar field to understand the underlying mechanism. For an arbitrary moment of 
time, the mode function for the scalar field is expanded as 

"fwAt > 0) = E <lnlnM \l r ^\ sin(n.7rg) sin(^y) sin(^z) , (55) 

V ^2x — ^Ix 
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where q{t) — — T~- If substitute it into the wave equation and use the orthogonahty relations, we get 



where Wm^„ n.(0 — ( "itt-xH-^^x — Li^) + Uy/ Ly + n^/ ) and the new antisymmetric coefficient is given by 

mInl(l + (-l)™==+"-) ^ 

gn^m^ = ^ 5 lor irix ^ rix- 

— 

There are two differences between pSI) and (|56p . Firstly, the antisymmetric coefficient vanishes when + is an 
even number when the cavity oscillates as a whole (given below ([35])). However, it vanishes when rrix + is an odd 
number when the cavity oscillates symmetrically. Secondly, the term „ „ in the left hand side of the former one 
is constant while it is time dependent for the latter one. The time dependent character of it gives a modification of 
the equation ([36]) 



^ nsci^rrix 

-2n'sm{nt) ^ gnxm^attlif^ (57) 
The second term in the right hand side is new. Let us study the the parametric resonance case, — ^ujuixnjjnz • 

Then, 

MSA gives equations for A{t) and B{t) 

7 ? 

Q Akxkyk^ , _2___x_f>kxkyk^ i Ak^kyk^ _ n /ro\ 

^-r^m^nyTi, ^ j-2 , , ^rnxfiyn, ^ (nxnyriAArrixnynA -^rixnyn, ~ ^ W7 
^x ^Tn^ 

2 2 

— r) fi^!c^y^' Akxkyk^ — nkxkyk^ _ q 

'-'T-i-'nixriyn^ jj2 ^ ^m^rLyn^ (nxnynz),('"x"»"i) "-xriyn^ " V^'^ ) 

where GV , was defined in (14111 and rix is a positive integer which satisfy the following relation 

^nxriyii^ — '^^nixriyn^ ■ (60) 

Let us give some examples. Firstly, we are interested in the uncoupled modes. For example, consider the mode 
(1, 1,0). If we solve (|58l59p endowed with the initial conditions ([55)1 . we get 

A}}°(r) = cosh(Ar) (61) 

Bll°ir) = -^=^ sinh(Ar) , (62) 



where A — —j= — . The number of photons is calculated (|24p 

< iVi,i,o >= sinh^(Ar) . (63) 

Compare this result with the one obtained in [18] where the authors assumed that only one wall is oscillating in the 
parametric resonance case. They found that < A^i,i,o >= sinh^(ADT), where Xd = A/2. 

Number of generated photons are the same for the following two systems: The cavity with single oscillating mirror 
with 2e and the cavity with two oppositely oscillating mirrors with e. Instead of increasing e by factor 2, the static 
wall is allowed to oscillate as described above. 

As a second example consider the mode (1, 1, 1). In this case, according to ([SO)) , the mode (5, 1, 1) is coupled. The 
solutions of (|58l59p become 

Alllir) = -^=S,,k^cosHX^T) (64) 

V ^^511 

Blllir) = -^^<55,,^sinh(Air), (65) 
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(t) = 0.681^1^ (sinh(Air) - sinh(A2r)) + -j^^ cosh(A2T) 
^''fiY (^) = 0.681-J^(cosh(A2r) - cosh(AiT)) - /"'^ sinh(A2T) (66) 

TT 25 TT 1 

where Ai — ^= and A2 ~ Then, the number of photons 

(Nixi) = sinh2(A2T) +0.1549 (cosh2(A2r) + cosh2(Air) - 2cosh2(A2T) cosh2(AiT)) 

(7V5,i,i) - sinh2(AiT) . (67) 

As a result, compared to the result obtained for a single mirror, the radiated photon flux is enhanced. 



VI. DISCUSSION 

Kim, Brownell and Onofrio proposed an experiment for the detection of the dynamical Casimir effect [isj . They 
considered a three dimensional cavity with a single moving boundary. Taking into account the limitation by the 
photon leakage of the cavity expressed through the optical quality factor Qopt, which saturates at the hold time 
T = Qoptl^i they gave formula for the maximum photon population for the parametric case and uncoupled mode as 

< iV>==sinh2 (2 g^pt ei) , (68) 

We have shown that if the two walls are moving opposite to each other with the same frequencies and amplitudes, e 
should be replaced by 2e. Hence 

< iV sinh^ (4 g^pt e t) , (69) 

The number of generated photons in the cavity is very sensitive to the product e Qopt- In current technology, the 
maximum value e Qopt = 1 The equation (|68p gives the number of generated photons A'^ = 13 if only one wall is 
in motion. If two walls move symmetrically with respect to the center of the cavity, (|69p gives the number of generated 
photons N = 745. 

The difference between the two cases is great if e Qopt — 2, which may be possible in the future. In this case, the first 
formula (l68|) gives N — 745. However, the second one ((69|) gives large number of photons N — 2.2 10®. 
From the experimental point of view, we think that the systems with two moving walls will play an important role 
for the detection of generated photons. 
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